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Abstract. We study the deformation caused by the influence of extrinsic curvature on a 
vacuum spherically symmetric metric embedded in a five-dimensional bulk. In this sense, we 
investigate the produced black-holes and derive general characteristics such as their masses, 
horizons, singularities and thermal properties. As a test, we also study the bending of light 
near such black-holes analyzing the movement of a test particle and the modification caused 
by extrinsic curvature on its movement. Accordingly, using the asymptotically conformal flat 
condition for the extrinsic curvature, an analytical expansion of a set of n-scalar fields can 
be defined and we show that the corresponding black holes must be large and constrained in 
the range of allowed values —1/2 < n < 1.8. As a result, they are locally thermodynamically 
stable, but not globally preferred. 
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1 Introduction 

In a previous publication [1], we studied the influence of extrinsic curvature on the rotation 
curves of galaxies obtaining a good agreement with the observed velocity rotation curves of 
smooth hybrid alpha-HI measurements. In doing so, we used the smooth deformation concept 
[1-4] based on Nash theorem of embedding geometries [5]. This allowed us the possibility to 
develop a more general approach to embedded space-times. Instead of assuming the string 
inspired embedding of a three-dimensional hypersurface generating a four-dimensional embed¬ 
ded volume, we look at the conditions for the existence of the embedding of the space-time. 
The embedding of a manifold into another is a non-trivial problem and the resulting embedded 
geometry resulting from of an evolving three-surface must comply the Gauss-Codazzi-Ricci 
equations. For instance, many of the currently brane-world models fail to comply with those 
equations. 

In the Randall-Sundrum brane-world model (RS) [6, 7], the space-time is embedded 
in a five dimensional anti-deSitter space AdS5. It is assumed that the space-time acts as a 
boundary for the higher-dimensional gravitational field. This condition has the effect that 
extrinsic curvature becomes an algebraic function of the energy-momentum tensor of matter 
confined to the four-dimensional embedded space-time. In more than five dimensions that 
boundary condition does not make sense because the extrinsic curvature acquire an internal 
index while the energy-momentum tensor does not have this degree of freedom. To obtain 
the so-called Israel condition it becomes necessary to consider the brane-world as a boundary 
separating two regions of the bulk, and find the difference between the expression calculated in 
both sides. The Israel condition only follows after the additional condition that the boundary 
brane-world acts as a mirror: Only in this case the tangent components cancel and the 
components depending on extrinsic curvature remain. This can be seen in the original paper 
by W. Israel [8], or in the appendix B of reference [2]. 

We should add that the Israel condition is not used in many other brane-world models, 
such as, for instance, the original paper by the Arkani-Hamed, G. Dvali and S. Dimopolous 
[9], sometimes referred as ADD model, and also the Dvali-Gabedadzi-Porrati model [10]. 
Rather, it was shown that the Israel condition is not unique, so different conditions may 
lead to different physical results [11]. Other brane-world models have been developed with 
no need of particular junction conditions, see [1-4, 12-15] and references therein. Moreover, 
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we should also add that such condition is particular to five-dimensional bulk spaces and is 
not consistent with higher dimensional embeddings. In the case of a Schwarzschild solution 
[20, 21], for instance, the embedding would be compromised in RS scheme and evinces the 
necessity of a more general framework. Unless we are restricted to a class particular models 
the Israel condition does not apply because only the mirror symmetry on the brane-world is 
dropped. 

In principle general relativity is exempt of the discussed Riemann curvature ambiguity 
problem because the ground state of the gravitational field, the Minkowski space-time, is well 
defined by the existence of Poincare symmetry (in special the translations). Such situation 
between particle physics and Einstein’s gravity was shaken by the observations of a small 
cosmological constant, then the emergence of the cosmological constant problem. This leads 
us back to the necessity of the embedding of space-times as the only known solution of the 
Riemann curvature ambiguity. 

The general solution of the embedding problem (mainly on Gauss-Codazzi-Ricci equa¬ 
tions) is shown in the paper by J. Nash on differentiable embedding [5], and generalized by 
Greene [22] including positive and negative signatures. The application of Nash theorem to 
gravitation is a landmark of our paper, providing a new tool to the gravitational perturbation 
issue. The embedding mechanics presented here is a legitimate condition in all other more 
general settings. 

The aim of our present work is to focus on the classical thermal stability of black holes 
embedded in five-dimensional bulk space-time. Recently, the study of the high dimensional 
space-times and the prediction of a new black objects (black holes, black strings, black rings, 
and so on) have been the focus of very active research [16-19]. Specihcally, it is worth noting 
that a Schwarzschild embedding was a problem discussed a long ago [20, 21, 26-28] and it is 
known that this geometry is completely embedded in six-dimensions [29]. In particular, the 
effects of a five dimensional bulk space-time endowed with an embedded spherically symmetric 
metric is still an open research arena since we do not have a closed analytical solution [30-33] . 
Notwithstanding, it may lead to new physical consequences in a search of a full description 
and understanding of astrophysical objects, mainly on stellar dynamics and gravitational 
collapse. Moreover, in this work we explore the consequences of such embedding by using the 
asymptotic condition for extrinsic curvature in order to attenuate such restrictions. 

The paper is organized as follows. In section 2, we briefly present a summary of the 
theoretical background of this paper restricted to the five dimensional space-time. In section 3, 
we study a spherically symmetric metric embedded in five dimensions. Moreover, in section 4, 
we analyze the constrains on the solution and the resulting horizons, which will be particular 
important to the thermodynamic analysis as presented in section 5. Finally, in the conclusion 
section, the final remarks and future prospects are presented. 


2 The induced four dimensional equations 

The present paper is not M-theory/Strings related. As well-known, the Randall-Sundrum 
(RS) model [6, 7] is a brane-world theory originated from M-theory in connection with the 
derivation of the Horava-Witten heterotic E8 x E8 string theory in the space AdS5 x S'5. In 
that framework, one obtains a compactihcation of one extra dimension on the or bifold 51/^2 
by using the Z2 symmetry on the circle 51. The presence of the hve dimensional anti-deSitter 
AdS5 space is mainly motivated by the prospects of the AdS/GFT correspondence between 
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the superconformal Yang-Mills theory in fonr dimensions and the anti-deSitter gravity in five 
dimensions. 

In a different approach, our proposal is initially based on the possibly to explore em¬ 
bedded space-times in a search of obtaining a more general theory based on embeddings (the 
brane-world models, of course, can be an example). Instead of assuming the string inspired 
embedding of a three-dimensional hypersurface generating a four-dimensional embedded vol¬ 
ume, we look at the conditions for the embedding of the space-time itself. 

Nash’s original theorem used a flat D-dimensional Euclidean space but this was soon 
generalized to any Riemannian manifold, including those with non-positive signatures [ 22 ]. 
For simplicity, we refer the term “Nash theorem” as valid also for pseudo-Riemannian man¬ 
ifolds. Another important aspect is that the seminal work on perturbation of geometry was 
proposed posthumously in 1926 by Campbell [23] using analytic embedding functions and 
later extended to non-positive signatures [24]. Although the Nash theorem could also be 
generalized to include perturbations on arbitrary directions in the bulk, it would make its 
interpretations more difficult, so that we retain Nash’s choice of independent orthogonal per¬ 
turbations. It should be noted that the smoothness of the embedding is a primary concern 
of Nash theorem. In this respect, the natural choice for the bulk is that its metric satisfy the 
Einstein-Hilbert principle. Indeed, that principle represents a statement on the smoothness of 
the embedding space (the variation of the Ricci scalar is the minimum possible). Admitting 
that the perturbations are smooth (differentiable), then the embedded geometry will be also 
differentiable. 

As it happens, Nash showed that any embedded unperturbed metric Qfj^u can be generated 
by a continuous sequence of small metric perturbations of a given geometry with a metric of 
the embedded space-time 

dfiu — dfiu 4 “ ^iiua 4 “ ^y 9 ^ k^pa^vab 4 “ ' ' ' . (^' 1 ) 


or, equivalently 


_ 1 dgpu 

2 dy 


( 2 . 2 ) 


where is the initially unperturbed extrinsic curvature and y is the coordinate related to 
extra dimensions. Since Nash’s smooth deformations are applied to the embedding process, 
the coordinate y, usually noticed in ridig embedded models, e.g [ 6 , 7], can be omitted in the 
process for perturbing an element line (depending on what one wants to do there are many 
different ways to embed a manifold into another classified as local, global, isometric, confor¬ 
mal or more generally defined by a collineation, rigid, deformable, analytic or differentiable 
embedding). Differently from RS models, using Nash embedding framework, we can develop a 
model that considers the extrinsic curvature as a dynamical (physical) quantity [1-4, 12-15]. 
Also, the paper should not be confused with a paper on quantization deformation, since our 
arguments are all classical. 

To start with, we consider an example with a Riemannian manifold 14 with metric g^i,, 
and its local isometric embedding in a five-dimensional Riemannian manifold V 5 , given by a 
differentiable and regular map A : l 4 —t V 5 satisfying the embedding equations 


= g^iu, X^^pTi^Gab = 0, ri^rj^GAB = 1, A, B = I..D (2.3) 


where we have denoted by Gab the metric components of V 5 in arbitrary coordinates, and 
where fj denotes the unit vector field orthogonal to l4. Concerning notation, capital Latin 
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indices run from 1 to 5. Small case Latin indices refer to the only one extra dimension 
considered. All Greek indices refer to the embedded space-time counting from 1 to 4. In 
general we have D = n + 1 with the index D representing the bulk space dimension and the 
index n represents the embedded space dimension. 

The extrinsic curvature of Vn is by definition the projection of the variation of r] on the 
tangent plane : 

kfj.u = ^j^Qab = X^ Q AB (2.4) 

The integration of the system of equations Eq.(2.3) gives the required embedding map X. 
Moreover, one can construct the one-parameter group of diffeomorphisms defined by the map 
hy{p) : Vd —t Vd, describing a continuous curve a{y) = hy{p), passing through the point 
p G Vn, with unit normal vector a'[p) = r]{p). Thus, the group is characterized by the 

def def 

composition hy o h±yi (p) = hy±y' {p ), ho (p) = p. Accordingly, applying this diffeomorphism 
to all points of a small neighborhood of p, we obtain a congruence of curves (or orbits) 
orthogonal to Vn- It does not matter if the parameter y is time-like or not, nor if it is positive 
or negative. 

If one defines a geometric object ui in Vn, its Lie transport along the flow for a small 
distance 5y is given by 14 = -|- 5y£y(l, where £rj denotes the Lie derivative with respect to 

y. In particular, the Lie transport of the Gaussian frame {X^,fj^}, defined on Vn gives 


= X^,y + 5y £yX^,y = X^,y + 6y (2.5) 

7]^ = + 5y = y^ (2.6) 

However, from Eq.(2.4) we note that in general 7 ^ y^y. 

In order to describe another manifold, the set of coordinates need to satisfy the 
embedding equations similar to Eq.(2.3) as 

^^Qab = g^lv, Z^^yy^QAB = 0, y'^y^GAB = 1 (2.7) 

Replacing Eq.(2.5) and Eq.(2.6) in Eq.(2.7) and using the definition from Eq.(2.4), we obtain 
the metric and extrinsic curvature of the new manifold 

g^JLu = g^lu - “^ykyu + y^g^^'kypkua ( 2 . 8 ) 

kyy = kyu “ 2ygP'^kypkyn (2.9) 


Taking the derivative of Eq.(2.8) with respect to y we obtain Nash’s deformation condition 
as shown in Eq.(2.2). Moreover, the integrability conditions for these equations are given by 
the non-trivial components of the Riemann tensor of the embedding space^. 

^To avoid confusion with the four dimensional Riemann tensor Ra/s-yS, expressed in the frame {Z^,ri^} as 

^'R-ABCdZ^ ^o,Z^ A ~ Rap-yS +{ka'ikps — kaskpj) ( 2 . 10 ) 

^'R-ABCdZ^ ^aZ^ = ka[p,A ( 2 . 11 ) 

These are the mentioned Gauss-Codazzi equations (the third equation -the Ricci equation- does not appear in 
the case of just one extra dimension). The hrst of these equation (Gauss) shows that the Riemann curvature 
of the embedding space acts as a reference for the Riemann curvature of the embedded space-time. The 
second equation (Godazzi) complements this interpretation, stating that projection of the Riemann tensor of 
the embedding space along the normal direction is given by the tangent variation of the extrinsic curvature. 
Notice from Eq.(2.10) that the local shape of an embedded Riemannian manifold is determined not only 
by its Riemann tensor but also by Its extrinsic curvature, completing the proof of the Schlaefli embedding 
conjecture by use of Nash’s deformation condition in Eq.(2.2). This guarantees to reconstruct the five¬ 
dimensional geometry and understand its properties from the dynamics, in this case, of the four-dimensional 
embedded space-time. 
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As in Kaluza-Klein and in the brane-world theories, the embedding space V 5 has a metric 


geometry defined by the higher-dimensional Einstein’s equations for the bulk in arbitrary 
coordinates 


'TIab - -^T^Qab = , 


( 2 . 12 ) 


where the metric of the bulk is denoted by Gab and we have dispensed the bulk with a 
cosmological constant, since for the present application in astrophysical scale the induced 
four-dimensional cosmological constant has a very small value and can be neglected. The 
energy-momentum tensor for the bulk of the known matter and gauge fields is denoted by 
'^AB- constant a* determines, in this case, the five-dimensional energy scale . 

In what concerns the confinement, the four-dimensionality of the space-time is an ex¬ 
perimentally established fact, associated with the Poincare invariance of Maxwell’s equations 
and their dualities, later extended to all gauge fields. Therefore, all matter which interacts 
with these gauge fields must be defined in the four-dimensional space-times. On the other 
hand, in spite of all efforts made so far, the gravitational interaction has failed to fit into 
a similar gauge scheme, so that the gravitational field does not necessarily have the same 
four-dimensional limitations and accesses the extra dimensions in accordance with Eq.(2.1), 
regardless the location of its sources. 

We assume that the four-dimensionality of gauge fields and ordinary matter applies to 
all perturbed space-times, so that it corresponds to a confinement condition. In order to 
recover Einstein’s gravity by reversing the embedding, the confinement of ordinary matter 
and gauge fields implies that the tangent components of in Eq.(2.12) must coincide 

with ( SnGTfj^u), where is the induced four-dimensional energy-momentum tensor of the 
confined sources. As it may have been already noted, we are essentially reproducing the 
brane-world program, with the difference that it is a general approach and has nothing to 
do with branes in string/M theory. Instead, all that we use here is Nash theorem together 
with the four-dimensionality of gauge fields, the Einstein-Hilbert principle for the bulk and 
a D-dimensional energy scale a*. The confinement is set simply as = SttGT^i,, where 

Tfj^u represents the energy momentum tensor of the confined matter, and T*^ = Tgg = 0. This 
can be understood in more general terms as a consequence of the duality operations of the 
Yang-Mills equations. It follows that D f\ F* is a three-form that must be isomorphic to 
the current 1-form j*. This condition depends only on the four-dimensionality of the space- 
time [34, 35] consistent with a plethora of experimental facts. Even that any gauge theory 
can be mathematically constructed in a higher dimensional space, we adopt the confinement 
as a condition, since the four-dimensionality of space-time will suffice in our case based on 
experimentally high-energy tests suggest [25]. 

With the above remarks, we may re-write the components of Eq.(2.12), using Eq.(2.2) 
and Eq.((2.7)) together with the previous confinement conditions [2] obtaining 



(2.13) 

(2.14) 


where the term results from the expression of TZab in Eq.( 2 . 12 ), involving the orthogonal 
and mixed components of the Christoffel symbols 


Q^u — 9 ^ ^fipkya kfj^yh 2 i^K h ) g^y , 
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where h? = is the (squared) mean curvature and is the (squared) Gauss 

curvature. This quantity is therefore entirely geometrical and it is conserved in the sense of 

= 0 . (2.16) 

Therefore we may derive observables effects associated with the extrinsic curvature capable to 
be seen by four-dimensional observers in space-times. A detailed derivation of these equations 
can be found in [2-4] and references therein. Hereafter, we use a system of unit such that 
c = G = 1. 


3 Induced spherically symmetric vacuum solution 


In the present study, we focus our work on obtaining an exterior vacuum spherical solution, 
i.e, = 0. To this end, we start with the general static spherically symmetric induced 
metric that can be described by the line element 


ds^ = B{r)dt^ — A{r)dr^ — r^dO^ — sin^ 9d(j>^ 


(3.1) 


where we denote the functions A{r) = A and B{r) = B. Thus, one can obtain the following 
components for the Ricci tensor: 


B" 


1 B' /A' B' 


1 A' 


2B AB\A^b) r a 
A' B' 

^99 = — 1 + TTw ( —r + “FT ) + yr 


2A 


A B 


B" 


= sin^ ORee 
IB' (A! B' 


Rtt — —FTT + T“r( “r + “Fr-r 


2A 4 A 


B 


IB' 


r A 


where we have ^ = A' e ^ = B'. 

From Eq.(2.13), the gravitational-tensor vacuum equations (with = 0) can be written 
in alternative form as ^ 

Rfiv T '^Qg^iu — Qfiu 

where we use the contraction Q = g^'^Q^p- 

The general solution of Codazzi equations Eq.(2.14) is given by 


(3.2) 


kf,u = f^ig^iu {no sum on g) , 

Taking the former equation and the dehnition of Qi_ip, one can write 


Qfii/ — ffj.dfJ.i' faffigfiu „ j fa ( f( 


(3.3) 



where 
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Consequently, we can write in terms of as 

Qtiu = Ufj,gfj,u {no sum on n). (3.4) 

Since is a conserved quantity, we can find four equations from = 0 that can 

be reduced to only two equations: 


/ (/l + /2)(/3-/4)=0 ; 
l(/3 + /4)(/l-/2)=0 . 


and result in the condition 

flh = /2/4 . (3.5) 

A straightforward consequence of the homogeneity of Eq. (2.14) and the condition Eq.(3.5), 
the individual arbitrariness of the functions can be reduced to a unique arbitrary function 
a that depends on the radial coordinate. Hence, the equation (3.3) turns out to be 

= a{r)g^y , (3.6) 

with a{r) = a. 

Thus, we can determine the extrinsic quantities 

Qfiu — 3 q ! g^n , ( 3 . 7 ) 

Q = Tr{Q^,) = 12a^ , (3.8) 

and substituting in the gravitational-tensor equation Eq.(2.13), one can obtain the following 
rr and tt components: 


B” IB' /A' B' 
TB “ 4(B VK ^ 
B" IB' /A' B'\ 

M 4^ 7 


^ ITI' 2 . ^ . 

J "774 =9« 

— - ^ = —90^{r)B . 
r A ^ ' 


(3.9) 

(3.10) 


From these equations, we have 

A' _ B' 

a( “ “(b ■ 

thus, AB = constant . 

in the same way as the calculation of the very known Schwarzschild vacuum solution, 
we impose the contour 

lim A(r) = lim B(r) = 1 , 

r^oo r^oo 

in order to approach the metric tensor to Minkowski tensor as r —?• oo, we have 


A{r) 


1 

B(r) ■ 


Furthermore, using the (0,9) component, one can obtain 


B(r) 


1-1 -1— / a^(r)mdr , 

r r 


(3.11) 
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and 


-1 


(3.12) 


A{r) = [B{r)]-^ 


1 + - + - [ a^{r)r^dr 
r r I 


where A; is a constant. 

Since spherically symmetric geometry is very constrained in a hve-dimensional embed¬ 
ding [29]. Based on the behavior of extrinsic curvature at inhnity, we set the asymptotically 
conformal flat condition on extrinsic curvature as 


lim kni, = lim a{r) lim g^u ■ (3.13) 

It is important to stress that the Nash theorem per se does not provide a dynamical set of 
equations for extrinsic curvature but shows how to relate extrinsic curvature to the metric 
through a smooth process of deformation making the local embedding regular and completely 
differentiable. As linir^ao 9fiu where is Minkowskian M 4 metric, the extrinsic 

curvature vanishes as it tends to inhnity, so the function a{r) must comply with this condition. 
Thus, we can infer that the function a{r) must be analytical at inhnity and to attend this 
constraint, we choose the simplest option 


a{r) = 

n=i 


\/~«o 


(3.14) 


where the sum is upon all scalar potentials and the indices i and s are real numbers. The index 
n may represent a set of scalar helds. The parameter ao has cosmological magnitude with 
the same units as the Hubble constant and its modulus is estimated as 0.677 [1] that will be 
used hereon. In order to keep the right dimension to Eq.(3.14), we have introduced a unitary 
parameter 7 * that has the inverse unit of Hubble constant [1] and dehnes the cosmological 
horizon in Eq.(3.17). From the geometrical point of view, Eq.(3.14) will not produce umbilical 
points on this spherically surface leading to a local additional bending in space-time. 

Using Eqs.(3.11) and (3.14), one can obtain an explicit form of the coefficient B{r) given 


by 


B(r) = 1 - + 1 ) _ 9«o ^2(i-n) 

r ^ 7*(2n — 3) 

n=i ^ ' 


(3.15) 


In terms on the correspondence principle with Einstein gravity, we set K(9ao -|- 1) = —AM, 
which remains valid even in the limit when a —)• 0 in order to obtain the asymptotically flat 
solution. In addition, the related potential is given by 


<h(r) 


2M ^ 9ao 2 (i-n) 
r 7 ^ 7 *( 3 - 2 n) 


(3.16) 


evincing how the initially spherical symmetric geometry is modihed by the influence of the ex¬ 
trinsic terms. Thus, we can write the line element in Eq.(3.1) modihed by extrinsic curvature 
as 

ds‘^ = 4>(r)(it^ — ] dr'^ — r^dO^ — sin^ OdcjP' , (3.17) 

<l>(r) 

for all n 7 ^: |. Hence, one can obtain physical singularities with the calculation of the induced 
four-dimensional Kretschmann scalar K and hnd 

_g(r, ao/'y*,M,n) _ 

r®(3 — 2n)^((r — 2M)(3 — 2n)r‘^^ + 9(ao/7*)^^^)® ' 
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K = 











r 

ao 

n 

Mass 

g{r, ao,M, n) 

K 

0 

0 

0 < n < 6 

any 

0 

undefined 

0 

0 

0 

any 

> 0 

diverges 

0 

0.677 

0 < n < 6 

any 

> 0 

diverges 

> 0 

0 

any 

0 

0 

0 

> 0 

0 

any 

> 0 

0 

> 0 

> 0 

0.677 

-10 < n < 10 

> 0 

> 0 

> 0 


Table 1. A summary of relevant values for the term g{r, a^, M, n) 
in determining different behaviors of Kretschmann scalar K. 


which preserves at first an intrinsic, physical and irremovable singularity at r = 0 and n = |. 

Another situations can be found when considering the term g{r,ao,M,n). This term 
represents a very large polynomial function with 384 terms that were omitted here in its 
explicit form. Bearing in mind that the value of the extrinsic parameter oq has magnitude 
0.677 or can be zero (if considering a vanishing extrinsic curvature) and analyzing the behavior 
of g{r, oo, M, n), one can summarize in the following table other singularities, a vanishing and 
positive K. 

As we can see, from the second to the fourth row in table (1) for any different values of 
n, the term g(r, ao, M, n) diverges leading to an undefined K. The fifth and sixth rows show 
an expected result since Riemannian curvature takes over when we have a vanishing extrinsic 
curvature. Interestingly, the last row shows the influence of the distortion cause by extrinsic 
curvature since even with a (non-forming) black hole with zero-mass [36], local tilde effects 
can be found once we do not have as a result the Minkowskian flat space-time. The range 
—10 < n < 10 was found once a larger range induces a very fast growing and decaying of the 
exponents but not observed in solar scales. In the section 4, we will study the possibility to 
get a more constrained range for those values of n. 

For the sake of completeness, we test the bending of light near the horizons in this 
geometry warped by Nash’s embedding. Using the line element for spherical symmetry in 
Eq.(3.1), one can study the geodesic motion given as the solution of the following equations 

A(-a*) = o, 

{r^9) — sin 9 cos 9^'^ = 0 , 

UT 

— sin^ 9 ^) = 0 
dr 

-AP -b A“ -b -b sin^ 9^'^ = 0 , 

(3.18) 


where r is the proper time and ^ = x. In addition, considering that the movement of a test 
particle takes place on a plane 9 = 7r/2, we perform a change of variables U = r~^ and find 
the following equation of trajectory 


d^U 

d^ 


+ AU + 




(3.19) 
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which yields 


^+^=3Mt/+|: ' “r-U 

Moreover, one can obtain the deflection angle as given by 

2 dr 


^ r 


n=i T 


nr-A 


11/2 


- TT, 


where the parameter b is the minimal distance from the source. Therefore, 


^ ^ AM ^ 


9ao/7*r^ ^'^dr 


-|3/2 ■ 


We note that 




(s)'-l 


1 3/2 


=^■’^ (3-2n)[(0"-l 

j-2-2«-2Fi( 3/2, 1 — n, 2 — n, r^/6^) 
' 2(1 -n) 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


where 2 F 1 is the gaussian hypergeometric function. Thus, one can And a generic expression 
for the shift Acf) as 

A4> = AcpGR + A(j)d , (3.24) 

where Ac^gr = ^ is the classical result provided from general relativity, and the second 
term A(f)d which is given by 


Acpd 


9a 

27 * 



2(l-n) 

- 5n + 3) 


2Fi( 3/2, 1-n, 2-n,r^/6^) , 


(3.25) 


which measures how strong is the deviation that depends on how intense is the deformation 
on the local space-time. For the sake of completeness, if we neglect ao, one obtains the same 
result obtained from general relativity. 


4 Constraints on the allowed values of the n-scalar fields 

In order to study the possible horizons by the element line in Eq.(3.17), we must constrain the 
values of the set of n-scalar fields. We start with analyzing the behavior of the gravitational 
potential on solar scales. In general, black holes are objects placed in galactic scales ruled 
essentially by newtonian potentials (~ ^) and its smooth deviations in decaying (~ ^) 

and growing not large as ~ r^. The increasing of gravitational potentials of order of ~ 
possibly can lead to a topological deformation of the space-time [37] or quadratically, and 
so on, should be quite worrisome and inconsistent with the solar system scale that turns the 
possible values of n dramatically constrained. Bearing this in mind, we analyze the allowed 
values of n, since we already have a first insight from the analysis of Kretschmann scalar in 
table (1). 

The general expression for the horizons of Eq.(3.17) can be found when we set gu = 0, 
so one can And 

n=i 
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which is a poiynomiai equation of order (3 — 2n) that gives a ciass of horizons to be dehned 
for aii aiiowed vaiues of n. Moreover, the corresponding mass is given by 



E 


j.2n-3 , ^ 

3 - 2n 


(4.2) 


In order to avoid overloaded notation, the unitary parameter 7 * is dividing the parameter ag 
hereon. Hence, we obtain interesting cases for the allowed values of n since the value n = 3j2 
is discarded due to the fact that it diverges Eq.(3.15). 

For instance, the Schwarzschild horizon can be straightforwardly obtained when we set 
ao = 0 , which means that extrinsic curvature vanishes and the space-time turns out to 
be asymptotically Minkowskian. For n = 0, one can recover the Schwarzschild-de Sitter¬ 
like solution and n = 2 mimics Reissner-Nordstrom solution. It is important to say that 
the Reissner-Nordstrom-Iike solution the term ag is interpreted as a tidal charge but not 
originated from a electromagnetic held, since it is a vacuum solution. Interestingly, in this 
particular case, differently from RS models, our tidal charge ao is originated from the extrinsic 
curvature and has a new physics associated to it, since it has a cosmological magnitude. The 
solution when n = 1 alone is a modihcation of the Schwarzschild solution by a term 9ao. 

For n > 3, the solutions produce a very low newtonian potentials (~ ^ and on). On the 
other hand, for n < — 1 the potentials grows very fast. Both limits are not verihed in solar 
system scales. Accordingly, we can set that the parameter n in summation must vary such 
as —l/2<n<3 hereon. 

An interesting case can be found when we set n = 1/2 where we have found the equation 

r — 2M — aor^ = 0 , 


with cosmological horizon Tc = and black hole horizon ri^h = ^ 

over, they exist when 1 — 8 aoM > 0 for M = and both horizons coalesce at Vc = r^h = 2ao‘ 
Interestingly, taking eq.(3.15), one can set n = where /3 is a parameter, and hnd the 
same power law found in [38], which was used for an alternative model to explanation of 
the dark matter problem. Another curious fact from the richness of this model is that when 
n = 5/2, which mimics an asymptotically metric function for a Bardeen black hole [39]. 


5 Thermodynamical stability 


In order to study the classical stability of this model, we focus on the determination of 
heat capacity and the resulting free energy. To this end, we calculate the related Hawking 
temperature, which depends on how curved is the space near a black hole. Therefore, one 
can hnd the local surface gravity as 


1 dgi4 

2 dr horizon) 


^ 0^2 (1 -^) 

— yo^n ~/ -r 

rh ° (3 - 2n) 


r 


l-2n 

h 


Consequently, the Hawking temperature is given by 



47rr2 


[rh + 9a 


2 3-2n] 
O^h \ 


(5.1) 
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Figure 1. The heat capacity is plotted as a function of the horizon with the allowed values of n. In 
the left panel is shown the behavior for the range in the allowed values —1/2 < n < 1.8 (solid line) 
and for —1/2 < n < 3 (dashed line). In the right panel, we compare the case —1/2 < n < 3 (solid 
line) to Schwarzschild solution (dashed line). 


where we denote the event horizon r^. For the sake of completeness, we point out that when 
ao 0 we recover Schwarzschild with event horizon rh = 2M and Hawking temperature 

7^ _ 1 

— SwM ■ 

In addition, we determine the entropy defined by Bekenstein-Hawking formulae as 


S 


r^h I QM 


T dr 


dr , 


and using the corresponding mass in Eq.(4.2), one can find 


S = Trrl, 


(5.2) 


which is interesting since it preserves the positivity of the entropic growth. 

From the classical expression of heat capacity, one can find the alternative form in terms 
of the mass of the black hole 


/dM\ _ fdM\ 


/ drh 

\dT 


and obtain the resulting formula 


C = 2 'Kr‘i 


(l + Oa, 


2 ^2(l-n 




9a§(l -2n)r2(i-»^) - 1 


(5.3) 


In order to analyze the thermal stability, we plot the resulting figures of heat capacity in 
terms of horizon as shown in Fig.(l). 

The positive values of heat capacity indicate a local stability of the system against 
perturbations. As the plots indicate in the left panel Fig.(l), the initial constraints with 
— 1/2 < n < 3, as shown in dashed line, lead to negative values for heat capacity, which 
means that we find an unstable black hole. Interestingly, a very different situation happens 
to a more tight range at —1/2 < n < 1.8, as indicated by the left panel in Fig.(l) in solid line. 
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Figure 2. The present figures show free energy as a function of the horizon and temperature for the 
allowed values of n. In the left panel is shown the behavior for the range with the allowed values 
— 1/2 < n < 1.8 (solid line) and for the values —1/2 < n < 3 (dashed line) with free energy as a 
function of the horizon. In the right panel we compare free energy to temperature and is shown the 
behavior for the ranges for the allowed values —l/2<n<1.8 (solid line) and —l/2<n<3 (dashed 
line). 
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Figure 3. The present figures show free energy as a function of temperature T for very small black 
holes. In left panel is shown the behavior for the range for the allowed values —1/2 < n < 3 and in 
right panel for the case —1/2 < n < 1.8. 


For small black holes we have an unstable solution. At the point of the horizon around at 0.5 
we find the extreme points of heat capacity indicating a phase transition. Accordingly, after 
that point, the larger is the horizon the more stable is the black hole against perturbations. 
In the right in Fig.(l), we compare the case —1/2 < n < 3 (solid line) to the standard 
Schwarzschild solution (dashed line). It results that the solution in the range —1/2 < n < 3 
is more unstable than the Schwarzschild one. 

Another important aspect is to study the global strength of stability which is given by 
the free energy. When it is related to temperature, it can be a valuable tool to study phase 
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transition regimes. The Helmholtz free energy is given by 


F = M -TS , 


and with Eqs.(4.2), (5.1) and (5.2), one can obtain 

The lesser values of the free energy indicate a global stability and the occurrence of Hawking- 
Page phase transition [40, 41]. We have checked two cases considering the free energy as a 
function of the horizons and latter as a function of temperature, as shown in the left and 
right panels in Fig. (2). 

As the plots indicate in the left panel in Fig. (2), the initial constraints with —1/2 < n < 
3, as shown in dashed line, lead to negative values for the free energy and then to unstable 
black holes as the horizon grows. A narrow window for transitory stability occurs around at 
0.5 < < 1, which is compatible with the left panel in Fig.(l). A similar pattern occurs 

when considering the temperature as shown in the right panel in Fig. (2) for a hxed radius 
R = 10. The transition phase occurs only at the origin with T = 0, so we have the extremal 
case. For hot space, we do not have phase transitions due to the fact that the free energy 
drops to negative values, which leads to several issues [15, 42, 43]. Conversely, we have an 
unstable solution as the heat capacity indicates. 

It is important to point out that in this case the higher is the free energy, the lower 
is the temperature leading to negative values, with can be discarded if we neglect unstable 
solutions. The same conclusion we find for very large black holes (with horizon ~ 10^^). 
For very small black holes (with ~ 10“^^), as shown in left panel in Fig. (3), the free-energy 
is positive and tends to a hxed value of T = 10. Thus, the system is not thermically stable 
and transition phases do not occur. The conclusion is that in that range, it does not provide 
stable small black holes. 

In the second case with —1/2 < n < 1.8, as shown in solid line in the left panel in 
Fig. (2), a phase transition occurs around zero free energy, not much different from the region 
presented in the left panel in Fig. (3). As we consider the temperature, as we go to T = 0 in 
right panel in Fig. (2), the free energy goes to a maximum at zero which indicates a thermal 
equilibrium. It means that for R = 10, and so on, we will not have a phase transition which 
is compatible with the heat capacity behavior. As a result, a large black hole will not have 
a phase transition as it grows if the temperature can be kept around zero, which dehnes an 
extremal black hole. Differently from the previous one, this case provides negative values 
for free energy which can give a more global stable situation, but with the price that the 
temperature must be negative. The same situation we have observed for very large black 
holes (with 1012). 

On the other hand, for very small black holes with a radius of the order of lO”^^ we have 
observed that they reach stability next to T = 5. At below radius of the order of 10“^ and 
so on, they stabilize at T = 10, and the free-energy continues to drop down. For a horizon of 
the order of 10 “ fj-gg energy is quite small of the order of 10“ which leads the conclusion 
that only small black holes can be stable at a certain range of temperature 0 < T < 10. As 
an overall conclusion, the best solution occurs at —1/2 < n < 1.8 predicting large black holes 
locally stable, but not globally preferred. 
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6 Remarks 


In this paper, we have discussed the local and global thermodynamic classical stability proper¬ 
ties for a class of static black holes embedded in a five-dimensional bulk space-time. Applying 
Nash embedding theorem to a static spherical symmetric metric, we have found a modifica¬ 
tion induced by the extrinsic curvature. Due to the fact that this metric has a restricted 
embedding in five-dimensional bulk space, we have dealt with the homogeneity of Codazzi 
equation in Eq.(2.14) using the asymptotic condition on extrinsic curvature based on the 
smooth principle of Nash theorem. Hence, the metric perturbation and the embedding lead 
naturally to a higher dimensional structure. 

We have shown that an analytical solution for embedded spherically static black holes 
is obtained. Analyzing the modified line element in Eq.(3.17) we have obtained the related 
horizons, which we have studied the behavior of a light ray near such horizons. As a result, 
we have found that the more extrinsic curvature is felted the more is the deformation on the 
space-time and the deviation is stronger than that one obtained from general relativity. 

Concerning the horizons, we have got a set of scalar fields initially constrained by the 
potentials produced in solar scale. Moreover, the determination of thermal quantities for 
such black holes, we have found a tight constraint for such set of scalar fields varying at 
— 1/2 < n < 1.8. As an overall conclusion, we have obtained a restricted model of stable 
modified black holes (initially spherically symmetric) as compared to RS models. In this 
range, we have obtained a local stable behavior for large black holes with a phase transition 
at rh ~ 0.5. 

We also have shown that the global stability is not a preferred state possibly influenced 
by the restricted embedding in five-dimensions, which may seem not suffice for an appropriate 
description of the phenomena. Accordingly, in order to obtain a more general situation, both 
global and local stability, an embedding in a more larger bulk seems an unavoidable situation. 
The entropy is that one expected vrr^ that grows positively. In RS model solutions for a 
Schwarzschild back-hole seems to be unstable [44]. In addition, since our analysis is classical, 
mini black holes are not allowed. 

Differently from RS model, we obtained a richer general set of static Black holes, e.g, 
Bardeen-like black holes and Reissner-Nordstrom-like black-holes. In this particular case, our 
tidal charge ao is originated from the extrinsic curvature and has a new physics associated to 
it [1-3], since it has a cosmological magnitude. As future prospects, a more general study of 
including rotation should be made concerning the canonical and grand canonical ensemble in 
a higher dimension, which are in due course. 
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